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ABSTRACT 


This repoi"t summarizes the results of using automated algorithms 
for removing aliases in NOSS-SCAT data. The algorithms used for alias 
removal consist of histogram analysis, local averaging and curve 
fitting. The histogram analysis is used to determine the degree of 
homogeneity of the wind field defined by the largest probability alias 
vector at each grid point. Depending on the degree of homogeneity, 
either a local averaging technique or a curve fitting technique is 
used to find a "preferred" wind direction at each grid location. The 
alias directions are compared with the preferred direction at each 
grid location and one of the multiple aliases is chosen as the true 
direction based on its closeness to the preferred direction and its 
probability. 

The alias reissval algorithms were applied to a number of simulated 
NOSS-SCAT data sets. The alias removal algorithms produ'zed correct wind 
vectors with an accuracy of over 955^. This compares with about 80 to 
85% of the grid points at which the^ alias with the largest probability 
is also thi correct solution. / 
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1.0 INTRODUCTION 

The major steps in the alias removal procedure are shown in Figure 1. 
The raw data is rear anged first into 10x1 0x26 arrays where each array 
contains all the data from a processing area covering 10 grid locations 
along-track and 10 grid locations across-track. The 10x10 processing 
area size was determined from consideration of an area small enough 
in size to provide some degrea of homogeneity and large enough to 
contain enough information about underlying streamline patterns. The 
streamline patterns were determined using the alias direction with the 
largest probability at each grid location. (Most of the processing 
was done on the array of largest probability aliases, i.e., the primary 
alias at each location.) 

The primary alias wind field over the processing area was used to 
compare histograms of wind directions, v/ind sprad and wind speed 
gradient. For each histogram, the mean, standard deviation and the 
entropy were also calculated. The standard deviation and entropy were 
used to determine the degree of homogeneity or uniformity of the wind 
field over a given area. If standard deviation and entropy were small, 
then the local average value of the wind direction was used to aid in 
alias removal. At locations where the local standard deviation. was 
large, then the average value of the primary alias directions over the 
entire processing area was used for alias removal. 

Higher values of standard deviation and entropy indicates the 
presence of a nonuniform wind field over the processing area produced 
by convergent set of streamlines, circulation patterns near lows, col 
point, or a frontal zone. In these cases a least square curve fitting 
algorithm is used to fit a streamline pattern to the primary alias wind 
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Major steps of the alias removal algorithm 
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field. The streamline pattern establishes a preferred direction at 
each point 1n the context of the directions at other points within the 
processing area. 

Given a context-derived preferred wind direction at a grid location, 
a choice Is then made as to which one of the multiple solutions Is the 
correct solution at a given location. This choice 1s made on the basis 
of the difference between an alias direction and the preferred direction, 
and the probability of the alias direction. The final output of the 
algorithm Is a unique wind direction for each grid location. The accu- 
:’acy of the algorithm Is evaluated by counting the number of locations 
at which the algorithm selects, as Its output, the alias direction 
closest to the true wind direction. 

The alias removal algorithm was used to process a simulated NOSS- 
SCAT data set (Rev 0825 among others). Histogram analysis Indicated that 
the left track contained slowly changing wind fields and the local aver- 
aging technique was used for alias removal. This algorithm selected the 
correct alias, I.e., the alias closest to the true direction, at 448 
grid points out of 449 points processed. The largest probability alias 
was the correct alias at 436 out of 449 points. 

The right track of the data set contained two lows and .a col point. 
These locations were easily Identified using the histogram analysis. 

The curve fitting approach was tested In the areas containing the lows. 
The alias removal algorithm selected the correct aliases at 133 points 
out of 143 points processed. The largest probability alias was the 
correct alias at 114 out of 143 points. 

Similar results were obtained for a number of other simulated NOSS- 
SCAT data sets. Details of the alias removal algorithm are given in the 
following sections of this report. 
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2.0 REFOBHATTING OF THE DATA SET 

The fii"st operation that was performed was the reformatting of the 
data. The data for the two tracks were separated and each track was 
broken up Into 10x10 blocks. There was an overlap of 5 rows between 
each of these blocks. The size of the measurement vector for each bln 
was reduced from 37 In the original tape to 26 In the reformatted 
version. The components of the measurement vector chosen were: 

1. column coordinate of the bin, 

2. row coordinate of the bin, 

3. average latitude of the measurements in the bin, 

4. average longitude of the measurements in the bin, 

5. average x coordinate of the measurements, 

6. average y coordinate of the measurements, 

7. wind speed from the original (true) wind field iata, 

8. wind direction measured from north from the original (true) 
wind field data, 

9. average wind speed (average of speed associated with each 
alias), 

10. number of aliases in the bin, 
for each of the four aliases: ‘ - ” • 

a. alias probability, 

b. estimated speed for this alias, 

c. estimated direction for this alias, 

d. magnitude of the residual. 
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3.0 HISTOGRAH <>NALYSES 

A simple statistical analysis of each block was peirfomied. The 
variables chosen were: 

(1) average wind speed, 

(2) wind direction associated with the largest alias probability, 

(3) average wind speed gradient. 

Each of the above variables was characterized by the following 
statistical descriptors: 

(a) mean, 

(b) standard deviation, 

(c) entropy, 

(d) histogram, 

(e) minimum, flt^ximum and their locations. 

A weighted histogram was used for wind direction, 1,e., the count In a 

cell, Instead of being Incremented by 1 Is Incremented by the probability. 

The formulae used to calculate mean and standard deviation of the wind 

directions are described In Appendix I. The gradient operator used 

on the average speed was a "Max Gradient" and is described In Appendix 

II. The average speed (and hence the gradient) varied from 0 m/s to 

24 m/s. . To construcju.ihfi.,histogram.this:^*.aijge was .d1vided;into If - - 

cells. The wind direction range (0“ to 360®) was also divided into 16 

16 , 

cells. As a consequence the entropy varied from 0 to ]C -1^100216 = 4. 

An entropy of 4 is indicative of an uniform distribution, and 0 is indi- 
cative of a histogram that consists of one nonempty cell. 

The true and aliased wind fields over the processing area are 
shown in Figures 2. a and 2.b, and a summary of the statistical para- 







SEA CHICKEN PLOT-WIND FLOW IN DIRECTION AWAY FROM CENTER 


8 




meters of data within 10x10 processing areas are shown in Tables l.a 
and l.b. Specific examples of aliased wind fields, the histograms and 
the values of mean, variance and entropy are given for two processing 
areas in Figures 3. a and 3.b. The examples show clearly that the 
histogram analysis indeed provides an indication of whether the wind 
field is uniform or nonuniform over the processing area. Entries in 
Table l.b, for example, show that the entropies of the wind speed 
gradient and wind direction are significantly larger in processing 
areas 4-8 which contain the lows. Entropies lower than 2 seem to 
indicate uniform wind fields, while higher values correspond to lows, 
col points, etc. This information is used in the next step of the 
alias removal algorithms. 
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Alia!.‘ed 
Wind Fiold 


0.,8 5.3. ,9.8 14.3 18.8 23.3 

Histogram of Wind Speed Gradient 
Mean *1.2 ni/s 

Standard Deviation = 0.5 in/s 
Entropy - U.8 


IU3 .28.8 146.3 ^13.8 281.5 348.8 

Histogram of Wind Direction With 
Largest Alias Probability 

Mean « 320.8° 

Standard Deviation = 24.2° 

Entropy - 0.9 


Figure 3. a. Example of histograms (left track). 




Histogram of Wind Speed 

Mean = 9.4 m/s 

Standard Deviation = 5.1 m/s 

Entropy - 3.4 


Histogram of Wind Direction With 
Largest Alias Probability 

Mean = 358,2® 

Standard Deviation = 88.3° 
Entropy =3.7 


Figure 3.b. 


Example of histograms (right track). 
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4.0 ALIAS REMOVAL BASED ON AVERAGES 

The following conclusions can be drawn from the values of the 
statistical parameters for the processing areas in the left track 
shown in Table l.a: 

(a) There is not much change in wind direction from one block to 
another. The average wind direction for the whole track is 
3 4.05® with a standard deviation of 28®. 

(b) The entropy associated with direction is also small, (1.28 
on an average) indicating a 'spiky' rather than an uniform 
probability distribution.^ 

(c) The average wind speed gradient is small with an entropy 
close to zero. So the rate at which the average speed 
changes is low. This coupled with the fact that the wind 
direction varies slowly with a relatively small standard 
deviation (28®) is indicative of the absence of cyclonic, 
anticyclonic or other such patterns. 

Because of the gradual changes in the vHnd field along the left 
track it would appear that simple local averaging could be used for 
alias removal. So, a 3x3 window was chosen and the averaging was done 
^ver the 9.ientries tft.-±hfi..w.indow. .(Sinca-.WBrare working with angles r; 
care should be taken to follow the procedure outlined in Appendix I.) 
The direction associated with the largest alias probability was the 
variable we chose to work with. The local average at each point was 
taken as the preferred ;*ind direction and the alias closest to the pre 
ferred direction was selet as the correct alias according to the 
following algorithm: 
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where: 


k • max 
1 


d.r 


i'i 



0 ) 


Allas direction selected 0|^, 

d^ « angular difference between the preferred direction (local 
average) and the ith alias direction 0^, 
r^ * residue corresponding to a^., 

S « a positive constant, 0 < 0 < 1. 


The figure of merit used to measure the performance of the alias 
removal algorithm was the comparison between the ‘number of agreements 
between the alias closest to true direction and the alias closest to 
average direction* on the one hand and 'the number of agreements 
between the alias closest to'true direction and the alias with the 
largest probability* on the other. 

Local averaging by itself did not clean up the track satisfacto- 
rily, because there are cells where the direction v.'th the largest alias 
probability was not ve »7 close to the true direction and there were 
cells for which no measurement was available. This resulted in the 
wrong alias being chosen. And an error in one cell affects the average 
of the neighboring cells. 

To rectify this problem, for each Windbw not only is the local 
mean calculated but also the local standard deviation (cr), and if the 
standard deviation is above a certain value a then instead of the local 
average the average of the whole block (global average) is used as the 

average for that cell (see Figure 4a). This experiment was carried out 

for a - 5®, 20°, 30°, 45°, 60°. The performance was the same for all 

values of a that are given above. The correct alias was selected at all 

but one point in the left track: 
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Figure 4. a. Processing details. 


3x 3Window 

\ 10 X 10 Processing Area 
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Total number of cells processed « 449 

Number of times the primary alias 

was closest to the true direction » 436 out of 449 

Results of the alias removal 

algorithms * 448 out of 449 

The alias removal algorithm failed at one point, located at row 35, 
column 2, Detailed examination of Figure 4b shows that even though the 
alias removal algorithm failed at this poi'lt, the difference between 
the errors is small (8.7® compared to 31.3®). 



Figure 4.b. Alias directions at grid location 
(35,2), left track. 

Instead of a 9 cell average, an 8 cell average could also be 
used, i.e., for each 3x3 window omit the center cell from mean and a 
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calculations. Performance in this case was as good as in the 9 cell 
case. 

The averaging algorithm was applied to a number of other tracks and 
the results are shown in Figures 5, 6 and 7. In Figures 6 and 7, the 
location of errors in the primary aliases (with the largest probability) 
and the location of errors after alias removal are explicitly shown. The 
results shown in Figures 5*7 Indicate that the averaging algorithm indeed 
works very well, and the improvement in accuracy is of the order of 8-125S, 
from 36 % prior to dealiasing to about 98^ after alias removal. 
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Rev: 0825, Left Track 
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Figure 6. Hi las removal based on averages. 
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Figure 7. Allas removal based on averages. 
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5.0 ALIAS REMOVAL BASRD ON CURVE FITTING FOR STREAMLINES > RIGHT 

TRACK 

The wind field over- the right track area for Rev 0825 varies 
considerably. These variations are reflected in the values of the 
statistical parameters shown in Table l.b (also see Figure 3.b). For 
example the wind direction has a mean standard deviation of 56® com-* 
pared with 28® for the left track, and the entropy of tN wind speed 
gradient is 1,79 compared to 0.74 for the left track. These variations 
indicate that an average based alias removal algorithm will not be 
successful and the structure of the streamline pattern within each pro- 
cessing area has to be taken into account by the alias removal al- 
gorithm. We focused our attention on blocks 5, 7 and 8 in the right 
track and used a curve fitting algorithm to find the best fit stream- 
line pattern and use this pattern to aid in alias removal. 

But, before proceeding to the streamline pattern analysis method, 
the entire right track was first processed using a 3x3 window and an 
8 cell average. Next it was processed using a 9 cell average/global 
average with thresholding on a ® 30°. The results of these two opera- 
tions are summarized in Table 2 (for blocks 5, 7, 8). 

» There is a deteriaciLtion in performanceras we go from a simple 8. - 
cell average to the thresholding case. The results in Table 2 indicate 
that the global mean can not be used for alias removal since it carries 
no information on the structure of a cyclonic vortex. 

The curve fitting approach was considered next. In this approach, 
the streamline pattern in the vicinity of the low/cyclonic vortex was 
modeled by a set of concentric circles. This circular pattern was used to 




Block 

# 

Agreement Between Alias 
Closest to Average 
Direction and Alias 
Clijsest to True Direction 

Agreements Bc*-ween 
Alias Closest True 
Direction and < lias 
With Largest Probability 


8 Cell Avg. 

9 Cell/Global Avg. 
With a - 30® 


5 

35/47 

25/47 

33/47 

7 

38/46 

26/46 

36/46 

B 

45/50 

41/50 

45/50 


Table 2 

Results of Simple Local Averaging and Averaging 
With Threshold on Cyclonic Areas 

predict the wind direction at any point. In the cyclonic vortexes of 
blocks 5, '7, 8 the flow is toward tlie center. This was incorporated 
by adding a constant 'inflow angle' to the direction predicted by the 
circular curve fitting. These angle calculations are summarized in 
Appendix III. Experiments were performed with various inflow angles and 
the results will be described later. 

But at first in order to generate the circular pattern we need 
to locate the center of the cyclonic vortex. In order to do so we 
assume the center to be at a certain point. We then generate the wind 
directior4S at various points based on the circular pattern and find the 
error between these "circle fitted" directions and, say the true 
directions. The center is then moved to another point and the proce- 
dure is repeated. This generates a 3-dimensional map: location of 

center (row, column) versus rms error. Three such maps were generated: 

(a) maps of rms error between true wind directions and "circle 
fitted" directions, 

(b) alias direction closest to true direction and "circle fitted" 
directions. 
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(c) direction with largest alias probability and "circle fitted" 
directions (all with an Inflow angle of 15"). 

These three maps are shown In Tables 3.1, 3.2, 3.3 and 3.4. 

The cell with the minimum rms error gives us the location of the 
center. Examining Tables 3.1, 3.2 we see that rms error increases as 
we go away from the center as would be expected. In Table 3.3, this 
behavior, though present Is a bit erratic. Also, the locations of the 
center obtained from the three maps are different but adjacent to each 
other. In a real situation, only Table 3.3 would be available to the 
user. Using this table we can locate the cell with the minimum rms 
error and take this cell to be the center of the vortex. The streamline 
pattern corresponding to this center defines a unique, preferred wind 
direction at each grid location in the processing area. 

The alias removal algorithm Is now applied to each point In the 
processing area. A local mean and standard deviation are first cal- 
culated using 3x3 windows. If the local standard deviation Is less 
than 3 j°, then the local average Is taken to be the preferred direction. 
If not, the direction given by the streamline pattern is taken as the 
preferred direction. The alias closest to the preferred direction is 
selected as the estimate of the true wind direction. The alias proba- 
bility Is taken Into account 'In the selection process as described in 
the previous section. 

Several experiments were performed over blocks 5, 7, 8 In the 
right track to determine the sensitivity of the alias removal algorithm 
as a function of the location of the center of the low, and the inflow 
angle included, and the threshold on the local standard deviation. The 
performance of the alias removal algorithm was evaluated by counting 
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.COLUMN 


\ 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

48 

86.86 

1 78.52 

70.77 

64,52 

62.57 


68.13 

73.47 

78.3 

78.9 

47 

84.22 


6fs,42 

57.19 

53.43 


38.5 

63.91 

69.45 

73.03 

46 

82.43 


62.16 

52.0 

45.11 

44.92 

48.86 

55.48 

62.73 

69,35 

45 

82.31 

72.33 

61.45 

50.25 

40.63 

(MIN) 

37.84 

42,36 

50.36 

58.60 

66.49 

44 

83.17 

73.62 

63.35 

52.96 

43.38 

39,38 

43.89 

50,71 

58.44 

66,44 

43 

85.93 

77.62 

68.51 

59.23 

51.68 

48.3 

48.76 

54.28 

61.28 

68.08 


Table 3.1 

RMS Error Between "Circle Fitted" Directions 
and True Wind Directions (Block #5) 

(Inflow angle = 15°, a ■ 30°) 


\C0LUHN 

R0W\ 

18 

19 

20 

21 

22 

23 

24 

25 

26 

48 

86,7 


70.53 

64.65 

63.73 

63.89 

68.51 

74.81 

81.15 

82.95 

47 

83,79 


64.69 

56.9 

55.08 

55.09 

60.05 

66.86 

72.99 

77,32 

46 

82.04 


60.2 

50.12 

46.18 

46,15 

52.02 

59.42 

67.11 

72,83 

45 

81.53 

70.95 

59.07 

47.1 

39,96 

40.04 

45.75 

54.47 

63.03 

71.0 

44 

81.83 

71.69 

60.43 ' 

50.62 

40; 53 

(.MIN). 

39.02 

46,42 

54.6 

61. r 

69.77 

43 

83.67 

76.29 

65.48 

55.93 

48.29 

46.42 

50,83 

-1 

57.16 

62,79 

70.55 


Table 3.2 

RMS Error Between "Circle Fitted" Directions 
and Allas Closest to True Direction 

(Inflow angle = 15°, a = 30°) 
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\C01 

rom\ 

.UMN 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

48 

95.96 

89.10 

82.7 

77.46 

75.12 

7<1,5 

76.38 

77.05 

77.52 

79.75 

47 

93.49 

85.92 

78.29 

71.80 

68.56 

69.02 

72.03 

74.7 

75.68 

76.55 

46 

91,98 

83.40 

75,28 

69.08 

67.76 

67.24 

68.95 

71.5 

75.0 

76.85 

45 

91,34 

82,98 

74,58 

72.40 

71.54 

69.01 

67.77 

(MIN) 

68.79 

72.39 

76.92 

44 

91.63 

63,75 

77.03 

73.38 

71.66 

71.16 

65.86 

68,23 

71.09 

75.30 

43 

93.64 

86.75 

80.64 

75.65 

73.95 

U9.74 

66.09 

68.48 

70,97 

73,26 


Table 3.3 

RMS Error Between "Circle Fitted" Directions 
and Direction with Largest Allas Probability for Block #5 

(Inflow angle * a » 30") 


\C0L 

R0W\ 

UMN 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

60 

101.24 

92.3 

82.51 

72.37 

62.87 

55.57 

51.78 

48.77 

(MIN) 

51.39 

56.5 

59 

98.21 

88.95 

78.82 

68.39 

58.74 

54.47 

48.75 

46.73 

48.4 

53.32 

58 

96.05 

86.95 

77.01 

66.67 

61.9 

54.23 

50.38 

48.41 

48.53 

51,33 

57 

94.31 

85,78 

76.7 

67.93 

61,57 

58.28 

53.14 

47.4 

50.09 

50.26 

56 . 

9,2.77 

85.34 , 

,,77. 8 

.2P^.5,4„ 

66.03 

59 ..11- 

:54.71 . 

.50.01 

50.66 

52*56 


Table 3.4 

RMS Error Between "Circle Fitted" Directions 
and Direction with Largest Allas Probability for Block #8 

(Inflow angle = 0", o >= 30") 
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the number of times the algorithi.t selected the alias closest to the true 
wind direction. 

The results of the alias removal attempts on the right track are 
shown in Tables 4-7. These results, which correspond to g » 0 in 
Equation 1, indicate that the algorithms are fairly robust and show 
promise. As e was changed, the best results were obtained thus far 
with 3 = 0.5. In the './icinity of lows, the curve fitting approach 
selected the correct alias at 44 out of 47 points, and 89 out of 96 
points. While the overall accuracy of 133 out of 143 points may be 
improved further, it is considerably better than the accuracy of the 
largest probability alias which is correct at only 114 out of 143 points. 


Location 
of Center 
WJ.iii) 

RMS Error 
(From Table 
3.3) 

# of Agreement? 

Between Alias 
Closest to 'Preferred' 
Direction and Alias 
Closest to True 
Direction 

# of Agreements 
Between Alias 
Closest to True 
Direction and Alias 
Direction With 
Largest Probability 

45,20 

72.4 

35/47 


46,22 

67.24 

37/47 


45,24 

68. 79 

39/47 

33/47 

44,24 

68.23 

40/47 

33/47 

45,21 

71.54 

40/47 

33/47 

45,22 

69.01 

41/47 

33/47 

45,23 

67.77 

42/47 

33/47 

44,23 

65.86 (MIN) 

, ,1 

43/47 - 

33/47- 


Table 4 

Effect of Changing Center of Cyclone For Block 
Algorithm Used: 9 cell average/''circle fit" 
threshold on o = 30® 
inflow angle = 0° 






omGiriAt pmt fs 
OF poor! ooALirr 



Location 
of Center 

RMS Error 
(From Table 
3.4) 

# of Agreements 
Between Alias 
Closest to 'Preferred' 
Direction and Alias 
Closest to True 
Direction 

# of Agreements 
Between Alias 
Closest to True 
Direction and Alias 
Direction With 
Laroest Probability 

59,24 

46.73 (MIN) 


81/96 

58,24 

48.41 


81/96 

58,22 

54.23 


81/96 


Table 5 

Effect of Changing Center of Cyclone for Blocks 7 and 8 
(Blocks 7 and 8 contain a single cyclonic vortex with 
the center somewhere in Block #8) 

Algorithm Used: 9 cell average/"circle fit" 
threshold on a » 30° 
inflow angle =^j0° 



Inflow 

Angle 

# of Agreements 
Between Alias Closest 
To True Direction and 
Alias Closest to 
'Preferred' Direction 

# of Agreements 
Between Alias Closest 
To True Direction and 
Alias Direction With 
Largest Probability 


0 

43/47 

33/47 

5 

5 

43/47 

33/47 

5 

15 

44/47 

33/47 

5 

20 

44/47 

33/47 

5 

25 

44/47 

33/47 

5 

45 

38/47 

33/47 . . 

"7&8 ' 

' 6 


-'■“-81/96 

7&8 

15° 

88/96 

81/96 


Table 6 

Effect of Changing Inflow Angle 
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Value of 
Threshold (=a) 

# of Agreements 
Between Alias Closest 
To True Direction and 
Alias Closest to 
•Preferred' Direction 

# of Agreements 
Between Alias Closest 
To True Direction and 
Alias Direction With 
Largest Probability 

60® 

41/47 


45® 

42/47 



43/47 

33/47 

20° 

43/47 

33/47 

15® 

43/47 

33/47 

5® 

43/47 

33/47 


Table 7 

Effect of Changing Threshold Level For Block #5 

Locatives of errors before and after the alias removal algorithm 
was applied are shown in Figures 8, 9, and 10. Figure 10 displays a 
situation in the vicinity of a col-point to which the curve fitting 
(Appendix IV) approach was used with a good degree of success. While the 
averaging technique improves the accuracy of picking the correct solution 
from 85% to 95% and upwards, the curve fitting approach improves the 
accuracy from about 77% to 93% in the vicinity of singular points. Con- 
sidering the large number of errors in the prime aliases near singular 
points, this improvement in accuracy is jei^^significant. ^ 






2Q 

Figure 8. Allas removal using curve fitting. 


OF POOR 


Before De-aliasing 


Number of cells processed: 47 

Number of errors in 

the primary alias : 14 (33/47) 

Number of errors 

after processing : 3 (44/47) 
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Figure 9. Alias removal using curve fitting. 


ORIGINAL PAGE IS 
OF POOR QUALITY 




Location of Errors 
Before De-aliasing 


Location of Errors 
After De-aliasing 


Number of cells processed: 96 

Number of errors in 

ttie primary alias : 15 (81/96) 

Number of errors 

after processing : 8 (88/96) 
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Figure 10. Allas iemoval based on curve fitting, near a col-point. 


-V ♦ • ® + ©© 

..L4^.-V+f 4®++^ 
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-^4^> -V -V t » ♦ ^ ^ ^ 

■V "V -V ■«■ ' s i \ 

*V*V'V + ^ V-H \ 


Errors in direction with 
largest alias probability 


Location of 

errors after processing 
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6.0 CONaUSJONS 0R19INAI PAGE IE 

.. u ...... Of POpn OUAUTY^ . . ^ 

We have presented In this report a set of algorithms that can be 

used for removing aliases In the multiple wind vector solutions produced 
from NOSS-type scatterometer data. The algorithms are based on histogram 
analysis, local averaging and curve fitting. The histogram analysis Is 
used to determine the degree of homogeneity of the wind field defined 
by the largest probability alias vector at each grid point. Depending 
on the degree of homogeneity, either a local averaging technique or a 
curve fitting technique Is used to find a "preferred" wind direction at 
each grid location. The alias directions are compared with the preferred 
direction at each grid location and one of the multiple aliases is chosen 
as the true direction based on Its'closeness to the preferred direction 
and its probability. 

The alias removal algorithms were applied to a number of simulated 
NOSS-SCAT data sets. The alias removal algorithms produced correct wind 
vectors with an accuracy of over 95 %. This compares with about 80 to 85% 
of the grid points at which the alias with the largest probability Is also 
the correct solution. 

Additional Improvement In performance may be realized by combining 
our procedure with other algorithms. To help Improve our performance 
near singular points and frontal zones, a catalog of histograms for 
various categories of streamline patterns should be developed. 






/ 
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Appendix I 

Mean and Standard Deviation of Angles 


(a) Compute Histogram 



^2 

^8 

^9 



•=16 


®2 

®8 




^16 


22.5® 45® 180® 

■ (weighted) count in the ith cell 
^ mid point of cell i; « 11.25® + (i-1) 22.5® 


(b) Compute two different means 


8 „ 8__ 

i»l i“l 

IG 16 

i=9 1=9 

8 

Let X, = 

1=1 ^ 
16 

^2=ES 

i=9 

(c) If ’ Qy < 180® 


(mean of upper quadrants) 


(mean of lower quadrants) 


6 = 


®u^l ®u^; 
+ X 


6^ - 0y > 180® 


0 = 


ByXi 4 (?^ - 360)X2 


Xi + Xg 


where ? is the required mean. 
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(d) Variance 



where * m 1 n{| 0 ^-i’l »36O-l0^-6l) 




*k 
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Appendix II 
Max 6rad’';ent 

J“1 J j+1 



^1,j * average wind speed at location 1,j. 

To find We max gradient at location 1,j the following four masks are 
used on the above v/indow; 



Max gradient at (i.j) = MAX(|0, |. lOjl, iDj] . |D^| ). 
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Appendix III 

Modelling the Cyclonic Vortex as Concentric Circles 



Figure III. a. 


Modelling the cyclonic vortex as concentric circles. 

(The true north is offset by a® from the perpendicular. ) 



Figure Ill.b. Angle calculations in the various quadrants. 
The center of the cyclone is assumed to be at (0,0) 



i{) is the required wind direction measured with respect to the 
true north. In Quadrant I: (j> = 180 - e - a 

II: (ji = 0 - a 

III: !|) = 360 - 0 - a 
IV: 4 . = 180 + 0 - a 
To incorporate the inflow angle: 
a ■<- a + inflow angle. 


Appendix IV 
Model For Co1 -Points 



Figure IV. a. Curve- fitting model for col-points. 


